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0\ ' ABSTRACT 

(N 

Context. We present an evolutionary disc model of the edge-on galaxy NGC 5907 based on a continuous star formation history and a continu- 
ous dynamical heating of the stellar subpopulations. 

Aims. This model explains the disparate two observational facts: 1) the exponential vertical disc structure in the optical and NIR of the non- 
obscured part of the stellar disc and 2) the FIR/submm luminosity enhanced by about a factor of four near the obscured mid-plane, which 
requires additional dust and also stellar light to heat the dust component. 
OO ■ Methods. We use multi-band photometry in U, B, V, R, and I- band combined with radiative transfer through a dust component to simultane- 
(~) ously fit the vertical surface-brightness and colour index profiles in all bands adopting a reasonable star formation history and dynamical heating 
\& . function. The vertical distribution of the stellar subpopulations are calculated self-consistently in dynamical equilibrium and the intrinsic stellar 
■ emissivity is calculated by stellar population synthesis. 
Results. The final disc model reproduces the surface-brightness profiles in all bands with a moderately declining star formation rate and a 
O . slowly starting heating function for young stars. The total dust mass is 5.7 x 10 7 Mq as required from the FIR/submm measurements. Without 
I ' a recent star burst we find in the midplane an excess of 5.2-, 4.0-, and 3.0-times more stellar light in the U-, B-, and V-band, respectively. The 
O , corresponding stellar mass-to-light ratios are 0.91 in V- and 1.0 in R-band. The central face-on optical depth in V-band is = 0.81 and the 
-4— » ■ radial scale length of the dust is 40% larger than that of the stellar disc. 

Conclusions. Evolutionary disc models are a powerful method to understand the vertical structure of edge-on galaxies. Insights into the star 
. . formation history and the dynamical evolution of stellar discs can be gained. FIR/submm observations are necessary to restrict the parameter 
^ ■ space for the models. 

>3 : 

5_j Key words. Stellar dynamics - Galaxies: kinematics and dynamics - Galaxies: interactions - dark matter 

_ _ i 

1 . Introduction dex distribution. Therefore detailed models are necessary to ex- 

tract the intrinsic stellar disc properties from multi-colour pho- 

1.1. Vertical structure of stellar discs tometry 

The vertical structure of the stellar disc in late-type galaxies T n a ser i e s of papers van der Kruit & Searle (fT98ll [l982al 
is mainly determined by the star formation history combined |1982bi started to investigate the vertical and radial properties of 
with the dynamical heating er(i) of the stellar population due me surface-brightness profiles of edge-on galaxies. They deter- 
to gravitational perturbations. Young stellar subpopulations are m i ne d for a number of nearby galaxies (including NGC 5907) 
confined tightly to the midplane with small velocity dispersion, the radial scale length, the radial cutoff, and the vertical scale 
whereas the older subpopulations are distributed over larger length. They recognised that the shape of the vertical surface- 
heights due to their higher velocity disp ersion This bas ic fea- brightness profiles can be well fitted by a sec/i 2 -model disre- 
ture is well-know n in the Galaxy (Wielenl 19771 Freeman[l99l] garding the central minimum in the case of strong dust lanes. In 
Edvardsson et al. 119931 e.g.). Dynamical heating leads to an rece nt decades much work has been done on the phenomeno- 
increasing mean age of the population with increasing height logic al ana lysis of edge-on galaxies (e.g. Guthrie lT992l Pohlen 
above the midplane. This results in an increasing mass-to-light e t a j 20021. 

ratio and vertical colour index gradients from blue to red. In 

• • , , • ■ j-j.t.ijuj'j.11. In our Galaxy the stellar subpopulations of the thin disc 

principle, these colour index gradients should be directly ob- 3 F F 

, , • j , . , . c . , j . are not homogeneously mixed. The vertical thickness of the 

servable in edge-on galaxies but unfortunately, strong dust ex- & 3 . . 

.... i ii j . . .-ii • subpopulations increases with age (Wielen & Fuchs 1988 1. 

tinction and reddening often dominates the vertical colour in- F F & v 1 ' 

Wainscoat et al. ( 1989 1 observed a similar vertical structure in 

Send offprint requests to: A. Just the edge-on galaxy IC 2531. For the analysis of the intrinsic 

Correspondence to: just@ari.uni-heidelberg.de physical structure of galactic discs we started to investigate the 
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vertical surface-brightness and colour index profiles of edge- 
on galaxies in order to analyse the age distribution of stellar 
discs (Wielen et al. 1992 Just et al. 1996). We used evolution- 
ary stellar population models for the intrinsic light distribution 
and include dust extinction in the radiative transfer for the ob- 
servable vertical surface-brightness profiles. We found that in 
edge-on galaxies the vertical colour index gradients above the 
innermost part, where extinction dominates, are a signature of 
the changing properties of the stellar population with height 
above the midplane. De Grijs and Peletier (2000 1 also found 
from a large sample of edge-on galaxies that the stellar popu- 
lation properties in the discs can be seen in the vertical colour 
index gradients above the extinction features. 

For NGC 5907 and other nearby edge-on galaxies Xilouris 
et al. ( 119971 IT9 99 1 used an extensive radiative transfer code in- 
cluding scattered light to model the two-dimensional light dis- 
tribution of the stellar discs. Using exponential intrinsic lumi- 
nosity distributions different for each band (B,V,I in the case of 
NGC 5907) and an exponential disc for the dust component, 
they determined radial and vertical scale lengths of the star 
light and the properties of the dust component. The resulting 
relative low dust mass corresponds to a gas-to-dust ratio of 810 
(corrected for the Helium contribution). In recent years FIR ob- 
servations offered a new, more direct way to deduce the dust 
properties of galaxies. Popescu et al. (2000 1 developed a so- 
phisticated temperature model of the dust component in edge- 
on galaxies in order to fit the FIR observations and applied it 
to NGC 891. With this method Misiriotis et al. J2001I deter- 
mined for NGC 5907 a dust mass 4 times larger than Xilouris 
et al. (IT9991 For the source of the FIR excess they adopted an 
additional young stellar component in a very thin disc, which 
is completely hidden by dust extinction in the optical and NIR 
bands. Recently Stevens et al. y()05) used a two-temperature 
dust model to fit the IRAS and SCUBA data for a sample of 
spiral galaxies. For NGC 5907 they determine a dust mass only 
three times larger than the Xilouris et al. value. 

1.2. Content of this paper 

In this paper we present an evolutionary stellar disc model for 
NGC 5907, which explains both observational results: 1) The 
exponential vertical profiles in all optical bands are reproduced 
by the outer stellar disc dominated by the older thin disc pop- 
ulation and are no longer just a set of independent fit functions 
in each band as in Xilouris et al. (I1999> . 2) The large amount of 
dust determined from FIR/submm observations by Misiriotis et 
al. {2001 1 exceeding the dust component found in Xilouris et 
al. by a factor of four is necessary in our model to obscure the 
additional light of the young stellar subpopulations with low 
mass-to-light ratio near the midplane of the disc. There is no 
additional recent star burst required to explain the additional 
star light needed for the dust heating. 

The basis of our disc model is a vertical cut with a self- 
consistent detailed distribution of stellar subpopulations. The 
vertical distribution of the stars is determined by the dynami- 
cal equilibrium of the stellar subpopulations described by the 
star formation history of the disc and by the dynamical heating 




Fig. 1. Deep R image of NGC 5907, sum of 4 exposures. The 
field of view is 15 x 13.25 arcmin. N is top, E is left, (cf Sect. [3} 



function (measured by the increasing velocity dispersion). The 
gravitational forces of the gas component and the dark matter 
halo are included. The intrinsic stellar light distribution is cal- 
culated from population synthesis models of simple stellar pop- 
ulations (SSP). The resulting vertical profile is then extended 
to an exponential disc up to the cutoff radius i? max by scaling 
the surface densities with constant scale heights. The intrinsic 
structure is not changed as a function of radius. We compute the 
observable vertical surface-brightness and colour index profiles 
by radiative transfer calculations with an exponential dust com- 
ponent, where the inclination of the disc with respect to the line 
of sight is fitted. 

The parameters of the basic vertical profile are chosen to 
match all constraints at some fixed radius Rq of the disc, where 
the disc is dominating the light and the bulge can be neglected. 
At that radius we fix the relative contribution of the gas and 
dark matter component Q g = S g /Stot and Qh = ^h/Stot 
to the total surface density S t ot up to a maximum height z max 
above the midplane. We will use Rq = lOkpc « 2R S with 
radial scale length R s of the stellar disc, which is near the 
maximum of Ed/Sh in the case of an isothermal halo. The 
scale height of the gas component relative to the effective scale 
height of the stars z g /zo and Sh = 0h/cre, the ratio of the ve- 
locity dispersion of the dark matter halo and the maximum ve- 
locity dispersion of the stars are also determined at i?o- 

In an iterative process the crucial fitting functions 
SFR(t a — t) and a(t) are optimised together with the global 
parameters: inclination i, radial scale length R s , effective scale 
height z s , dust distribution, and total stellar and dust mass. 
For comparison with observations we use a series of vertical 
surface-brightness (and colour index) profiles parallel to the 
minor axes in U,B,V,R,I bands, which we obtained from cor- 
responding deep photometric images of NGC 5907. The same 
type of intrinsic disc model was also successfully applied to 
the solar neighbourhood in order to analyse the connection 
between the local velocity distribution functions of the main 
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sequence stars from the Catalogue of Nearby Stars (CNS4, 
Jahreiss & Wielen fl997i and the vertical density profiles of 
these stars (see Just l200T1l2002ll2003t . 

In Sect. 12 we collect the basic parameters of NGC 5907, 
which are necessary for the basic scaling of our disc model 
to the observational data. In Sect. [3] we describe the building 
blocks of our disc model. Sect. |4] gives the iteration proce- 
dure with details on the scaling method, radiative transfer and 
the multi-colour analysis. In Sect. [5] our observations and the 
derivation of vertical colour index and surface-brightness pro- 
files are presented. Fig.^shows a deep R image (cf Section|5]as 
an example of for the high quality of data. In Sect. |6] we show 
the final comparison of the profiles, give the global parameters 
of our best model and discuss in detail the intrinsic structure of 
the stellar disc. A discussion of some crucial aspects follow in 
Sect.0 Sect.[8]contains the concluding remarks. 

2. Basic parameters of NGC 5907 

Here we discuss the basic scaling parameters of NGC 5907, 
namely distance, radial and vertical scale length and inclination 
of the stellar disc, and the masses of the stellar, gaseous and 
dark matter component. 

2. 1. Distance to NGC 5907 

There is some uncertainty in the distance determination of 
nearby galaxies from the Hubble expansion, because the pe- 
culiar motion of the galaxies is not negligible. In this work we 
use a distance of D = 11 Mpc with Vqsr = 817km/s from the 
RC3 catalogue (de Vaucouleurs et al. 1991 ) and a Hubble con- 
stant of Hq = 75 km/s/Mpc. The resulting scale conversion is 
1" = 53.3 pc and the distance modulus is Am = —30.2 mag. 

Another distance determination for NGC 5907 using the 
Tully-Fisher relation leads to a slightly larger distance of 
D=11.6 or 12.0Mpc fSofue fl996l rIW7l . In this method the 
main uncertainty arises from the conversion of the edge-on lu- 
minosity to a face-on luminosity. Zepf et al. (2000) used a dis- 
tance of D = 13.5 ± 2.1 Mpc from a combination of the R- 
band Tully-Fischer relation with a peculiar motion model and 
discussed also (due to the lack of resolved giant stars in the 
outer bulge region from HST observations) the possibility of a 
significantly larger distance. 

2.2. Disc parameters 

For the outer parts of the disc, where extinction is not sig- 
nificant, the brightness distribution of the disc of NGC 5907 
can be modelled by exponential profiles in the vertical and ra- 
dial direction. The radial scale length increases from 3.7 kpc 
in H-band (Barneby & Thronson 1994 1 to more than 5 kpc 
in the blue (van der Rruit & Searle 1982a 5.7 kpc in J-band, 
Xilouris et al. [19991 5.02 kpc in B-band). Van der Kruit & 
Searle also determined the cutoff radius of the exponential disc 
at i? max = 19.3 kpc, which has some influence on the inner 
parts of the radial surface-brightness profiles due to projection 
effects. The cutoff was not included in the model of Xilouris et 
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al. (TJ999 ) leading to an underestimation of the intrinsic scale 
length. 

The effective scale height of the disc is approximately in- 
dependent of radius, but depends strongly on the model used, 
i.e. additional components like thick disc, bulge and also on the 
radial profile for edge-on galaxies. From moderately deep pho- 
tometry van der Kruit & Searle found z s = 410 pc in J/band 
and Xilouris et al. z s = 340 pc in B-band, and Barneby & 
Thronson z s = 320 pc in H-band, which are relatively small 
values. With a One-disc-model for the very deep R-band pho- 
tometry Morrison et al. (1994) determined a scale height of 
z s = 467 pc. 

From the shape of the outer isophotes van der Kruit & 
Searle found an inclination of i = 87.0°, which is widely used 
as a standard value. Xilouris et al. determined independently 
a value of i — 87.2°. Morrison et al. d 1994b observed a stel- 
lar warp in the outer parts of the disc. Therefore the intrinsic 
inclination may be even higher, since the warp smears out the 
isophotes additively. 

Since in our disc model we assume only one stellar disc for 
all five observed bands, we redetermine the radial and vertical 
scale length, the cutoff radius and the inclination in the fitting 
procedure. 

2.3. Masses 

For the construction of the disc model we need the relative sur- 
face densities of the different components. The scaling to the 
absolute values is then a result of the surface-brightness pro- 
file fitting (cf Sect. 14.3b . and the total masses of the stellar and 
gaseous components follow from the radial extrapolation to ex- 
ponential discs. We use the total hydrogen mass Mh and con- 
vert it to the total gas mass M g which will be used for the deter- 
mination of the gravitational force and to calculate the gas-to- 
dust ratio F. For the dark matter halo we need only a reliable 
estimate of the local density in the disc. 

2.3.1 . Gas masses 

The atomic gas mass is determined from the 21cm luminosity. 
We use the standard value 

A/hi = 6.9 x 10 9 M Q (1) 

(Dumke et al. [19971 Stevens et al. 120051 . The molecular gas 
mass is much more uncertain. It is derived by converting CO 
measurements using the conversion factor X, which gives the 
corresponding Hi mass. Dumke et al. ( 1997 i determined indi- 
vidually for the galaxy NGC 5907 a reduced conversion factor 
X yielding M Ha = 0.9 x 10 9 M Q . The standard X-value would 
have led to 1.8 x 10 9 Mq, which demonstrates the intrinsic un- 
certainty of the molecular gas mass determination. We use the 
molecular gas mass redetermined by Stevens et al. ( 2005 1 of 

M H2 = 1.7 x 10 9 M Q (2) 

(rescaled to the distance of D — 11 Mpc). For the contribution 
of He to the total gas mass we add 40% to the total hydrogen 
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mass Mh = Mm + Mjj 2 ■ This leads to the total gas mass of 
M g = 1.4 M H = 1.2 x 1O 1O M with (3) 
M H = Mhi + M Ha = 8.6 x 1O 9 M . (4) 

The radial distribution of HI is flatter and more extended 
than the stellar light distribution and is not exponential. The CO 
distribution is much more confined to the inner parts than that 
of HI (Dumke et al. U997t . For simplicity we will also adopt for 
the scaling of the gas fraction an exponential disc, but allow for 
a scale length larger than the radial scale length of the stellar 
disc. We will use the same spatial distribution as for the dust 
component, which is determined by our fitting procedure. 

2.3.2. Stellar and dark matter masses 

The determination of the stellar disc and dark matter (DM) halo 
mass from the rotation curve is very uncertain due to the well- 
known disc-halo degeneracy. Sackett et al. ( 1994 1 constructed 
three-component mass models (HI, stellar disc, dark matter 
halo) to fit the Hi-rotation curve of Sancisi and van Albada 
(19871. They found equally good fits with disc mass-to-light 
ratios of 1, 2, and 4 in the R-band. The models of Sackett et al. 
with M/L = 1 and M/L = 2 yield a range of surface density 
ratios of stellar disc and halo of 

— —2... 5 with Rq = lOkpc, z max = 5.1 kpc. (5) 

Eh Ro 

Typical mass-to-light ratios of stellar discs in the V-band are of 
order unity and are dependent on the average age of the pop- 
ulation. In the solar neighbourhood we have M/Ly — 0.78 
and the extrapolation to the solar cylinder is M/Ly ~ 1.4 for 
the mass-to-light ratio of the surface density (from CNS4 data, 
Jahreiss & Wielen lT997l . Our final model yields M/Ly = 0.9 
and M/L-r, — 1.0. Therefore we will use the M/Zr = 1 
model of Sackett et al. (I1994> for the scaling of the DM halo. 

Sofue (1996 1 presented the joint rotation curve of CO 
and HI, which shows in HI a slightly higher and more pro- 
nounced maximum, which is harder to reconstruct with a low 
mass disc of M/Lr < 2. From H-band photometry Barneby 
and Thronson d 1 9941 investigate mass models including a flat- 
tened bulge and determined a small bulge with scaling radius 
0.23 kpc and total mass of M B = 9 x 10 9 M . The CO rota- 
tion curve of Sofue ( 1996 1 confirms the kink produced by the 
bulge at V C (R = lkpc) w 200km/s of model b) in Barneby 
and Thronson (I1994> . As a consequence the pollution of the 
disc luminosity with bulge light at radial distances larger than 
3 kpc is very small. Therefore we decided to neglect the bulge 
component in our analysis in order to keep the number of free 
parameters small. 

2.3.3. Dust mass and extinction 

Our aim is to construct a disc model with a dust component, 
which has a mass comparable to 



M d = 6 x 1O 7 M , 



(6) 



the value determined by Misiriotis et al. (2001 1. We use the 



was confirmed by Xilouris et al. (fT999l for NGC 5907 and 
which also goes into the conversion factor of extinction to dust 
mass. In our model we determine the spatial distribution of ex- 
tinction Ay(R, z) (see Eg. I24>. This will be converted to the 
dust mass distribution by 

p d = 0.175 x 1O~ 3 A V M pc~ 3 (7) 
with Ay in [mag/kpc]. This is the same conversion factor as 



used by Xilouris et al. (119991 and Misiriotis et al. (2001 1. 

3. Building blocks of the disc model 

The aim of this work is to construct a physical model of the 
stellar disc in order to reproduce the vertical surface-brightness 
and colour index profiles from U,B,V,R, and I-band observa- 
tions. Therefore we pay much attention to the internal structure 
of the stellar disc. The gas and dust component and the dark 
matter halo are modelled in a simple way. The bulge contribu- 
tion is neglected, because the number of free parameters and fit- 
ting functions would be approximately doubled with very small 
effect on the parameters of the disc. A thick disc component is 
excluded by very deep photometry and the faint stellar halo is 
below our observational limit (Morrison et al 1994 1. We ap- 
ply our model to those regions of the disc where the thin disc 
strongly dominates. 

In this section we describe the different ingredients nec- 
essary to construct the disc model and compute their intrinsic 
properties. The projection onto the sky and the fitting procedure 
is given in Sect.|4] 

3.1. Self-gravitating disc 

The backbone of the disc is a self-gravitating vertical disc pro- 
file including the gas component in the thin disc approximation. 
In this approximation the Poisson-Equation is one-dimensional 



d 2 $ s 



dz 2 



AirGp{z) 



(8) 



where p{z) is the self-gravitating density and $ sc if(z) is the 
corresponding potential. Since we want to construct the disc in 
dynamical equilibrium, the density of the sub -components will 
be given as a function of the total potential pj (<3>) and not of 
height z. In the case of a purely self-gravitating thin disc (with 
^scif = no external potential) the Poisson equation can 

be integrated leading to 



d$\ 2 /"* 



(9) 



Then the vertical distribution is given by the implicit function 
z($) via direct integration 







/ d$' 


8ttG / p($")d$" 


/o 


Jo 



-1/2 



(10) 



standard extinction law of Rieke & Lebofsky (1985 1, which 



If an external potential is included, an iteration process is nec- 
essary to solve for the vertical distribution. In order to avoid the 
iteration at that point, we model all gravitational components 
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by a thin disc approximation. We include in the total potential 
$ the stellar component $ s , the gas component $ g and the dark 
matter halo contribution $h 



$(z) = $ a (z) + $ s (z) + $ h (z) 



(11) 



In order to obtain the force of a spherical halo correctly in 
the thin disc approximation, we use a special approximation 
(see subsection l3.4t . The relative contribution of the stellar, the 
gaseous, and the DM-component to the surface density (up to 
\z\ — z mal ) are given by the input parameters Q s ,Q g ,Qh- 

3.2. Stellar disc 

We assume that the disc of NGC 5907 is, like the thin disc of 
the Milky Way, composed of a sequence of stellar subpopula- 
tions with increasing vertical scale height with age. This can 
be parametrised by the star formation history and the dynami- 
cal heating function. For convenience we use here the normal- 
isation to the central profile at R = 0. The stellar component 
is composed of a sequence of isothermal subpopulations char- 
acterised by the IMF, the chemical enrichment [Fe/H](t), the 
star formation history SFR(t a — t), and the dynamical evolu- 
tion described by the vertical velocity dispersion <j(t). Here t 
is the age of the subpopulation running back in time from the 
present time t a = 12 Gyr (which is the adopted age of the disc). 
We include mass loss due to stellar evolution and retain the 
stellar-dynamical mass fraction g(t) (stars + remnants) only. 
The mass lost by stellar winds, supernovae and planetary neb- 
ulae is mixed implicitly to the gas component. 

With the Jeans equation the vertical distribution of each 
isothermal subpopulation is given by 



p s ,i(z) = p s0 ,j exp 



-$(*) 



(12) 



where p s j is actually a 'density rate', the density per age bin. 
The connection to the SFR is given by the integral over z 



/oo 
-oo 



(13) 



The (half-)thickness h p (tj) is defined by the midplane density 
p s0 ,j through 

SFR(U-t 3 ) 



" M " ~ 2ftp(tj) 
The total stellar density is 

p s (z) = / p s ,jO)df , 



(14) 



(15) 



which then determines the potential $ s {z) via the Poisson Eq. 
(|8j. The stellar surface density E s is connected to the integrated 
star formation So by the effective stellar-dynamical fraction 

5off 



with 



p s (z)dz = g c sS 



(16) 



S =[sFRdt and 9eS = , 17 ) 



The effective scale height z s of the stellar disc is connected to 
the maximum velocity dispersion rr c of the subpopulations and 
the total surface density S tot via 



CyZa C 7 , 



(18) 



27rGStot 

where z c is the scale height of an isothermal component above a 
disc with total surface density Etot- The shape correction factor 
C z is of order unity and is determined at z = (3 ± 0.5) z c . 

The metalicity Z affects the stellar lifetimes, luminosities 
and colours of the subpopulations. Z shows a vertical gradient, 
since the mean age of the stellar population is correlated with 
the vertical distribution. In order to account for the systematic 
influence of the metal enrichment we adopt a moderate metal 
enrichment similar to the solar neighbourhood (Twarog 1980 
Edvardsson et al. 1993 1. We use the enrichment law for oxygen 

+yin(l + a— ) with Y = J^——^ (19) 



Z 



t a ) ln(l + a) 

with Z normalised to the solar abundance and the conversion 
law [Fe/H] = 2[0/H] as a simple analytic description to ac- 
count for the enrichment delay by SNla. This model is bor- 
rowed from a closed-box model with an n = 2 Schmidt-law 
for the star formation (Lvnden-Bell ll975l Just et al. U996i . We 
use a = 5, initial and present metalicity z s = 0.4, Z p = 1.13 
corresponding to [Fe/H] = —0.8 and [Fe/H] =0.1 with 
solar metalicity Zq = 0.018, respectively (see Fig.|2j. 

In the fitting procedure (Sect. 14. 1 i a pair of star formation 
history and heating function is selected to derive the intrinsic 
structure of the disc. The star formation history and heating 
function of the final model are shown in Fig.|2] 

3.3. Gas and dust component 

For the gas and dust component we use simple models to ac- 
count for the gravitation of the gas and the extinction of the 
dust. In the radial direction we adopt an exponential profile 
with scale length i?d = qdRs, where we allow for a differ- 
ence in the scale length of the stellar component by the factor 
<7d- For the extinction by the dust component we use a simple 
exponential profile with vertical scale height z&. 

The vertical profile of the gas component, which is used for 
the gravitational force of the gas, is constructed dynamically 
like the stellar component. The gas distribution is modeled by 
distributing the gas with a constant rate over the velocity dis- 
persion range <r(t) of the young stars up to a maximum age t g . 
By varying t g we force the scale height of the gas z g to the 
same value as that of the dust component z^- The surface den- 
sity of the gas is related to the stellar surface density by the 
ratio Q g /Q s — S g /E s , which is determined at the reference 
radius Rq in the fitting procedure (cf Sect.|4). 

3.4. Dark matter halo 

The halo does not fulfil the thin disc approximation. For a 
spherical halo we get the vertical component of the force to 
lowest order from 

d£h _ GMr£ 

dz ~ i? 2 R ' K ' 
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Fig. 2. The upper panel shows the normalised SFR/Sq as 
function of normalised age for the final model of NGC 5907. 
The dashed line is the constant SFR of a comparison model. 
The lower panel gives the corresponding heating function 
a(t) I a normalised to the final velocity dispersion a c (full line) 
and the chemical enrichment (dotted line). The dashed line is 
the heating function of the solar neighbourhood used for the 
comparison model. 

with r 2 — R 2 + z 2 and Mr is the enclosed halo mass inside 
radius R. Comparing this with the one-dimensional Poisson 
equation from the thin disc approximation (Eq. [S] integrated 
over z near the midplane to lowest order for small z) 



d$ 

— = AnGpoz 
dz 

we should use for the local halo density 



PhO 



Mr 
AttR 3 



(21) 



(22) 



which exactly corresponds to the singular isothermal sphere. 
Therefore we can use the thin disc approximation also for the 
halo, if we use the local halo density pho an d the halo velocity 
dispersion o"h estimated from the rotation curve by adopting an 
isothermal spherical halo. The effect of a cored halo, anisotropy 
and flattening is neglected here. For other halo profiles cor- 
rection factors would be necessary introducing some inconsis- 
tency in the halo profile description and leading to a different 



local halo density. The latter would be more important, because 
the effect of the halo potential on the disc is stronger than the 
adiabatic contraction of the halo in the disc potential. 

3.5. Radial structure 

The parameters of the self-consistent vertical profile described 
in the last section are determined at a scaling radius Rq = 
lOkpc. The determination of the disc properties at Rq are de- 
scribed in Sect. 14. II For the radial structure of the disc we use a 
simple exponential extension up to a cutoff radius i? max - We do 
not account for the radial variation of intrinsic disc properties, 
which are 1) the star formation history or metalicity to fit the ra- 
dial colour index gradients, 2) the increasing dark matter mass 
fraction with radius, 3) the possible variation of the star/gas 
surface density ratio due to different radial scale lengths, 4) the 
breakdown of the thin disc approximation in the innermost part 
of the disc, 5) the bulge potential and luminosity in the inner 
region. 

The stellar disc model is then given by 



p s (R, z) = p s (z) exp(-i?/.R s ) 



for R < R„ 



(23) 



with radial scale length R s , cutoff radius i? max and constant 
scale height z s . 

For the gas and dust component we allow for a different 
radial scale length R g = i?d = qdRs with scaling factor q^. 
For the conversion of total masses to density distributions we 
do not apply a cutoff radius. The corresponding distribution of 
extinction is 



Ay(R, z) = Ay,o exp 



( ~ R \ z \ 



(24) 



\q d R s Zd, 
with central extinction coefficient Ay,o- 

3.6. Stellar population synthesis 

The intrinsic luminosities in the different bands are determined 
by stellar luminosity synthesis. We use the stellar population 
synthesis code PEGASE (Fioc & Rocca-Volmerange ll997t to 
produce "pseudo" simple stellar populations (SSPs), which are 
used for computing the intrinsic luminosity distribution of the 
stellar disc. This means that the PEGASE code is used to cal- 
culate the integrated colour indices for a stellar population cre- 
ated in a single star-burst at different time-steps. These SSPs 
are then used to assemble a stellar population with a given star 
formation history, in the sense that the star formation history 
is assembled by a series of star-bursts. Since we do the assem- 
bling of SSPs ourselves, we have to take into account our own 
treatment of chemical enrichment. 

Our "pseudo" SSPs are modeled by a constant star forma- 
tion rate with a duration of 25 Myr. As a input parameter for 
the PEGASE code we use a fixed Scalo-like IMF (Scalo [T986t 
given by 



dN 



M- a dM 

1.25 0.08 < M/Mq < 1 
2.35 for 1 < M/Mq < 2 

3.0 2 < M/Mq < 100 



(25) 
(26) 
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Fig. 3. colour index evolution of single stellar populations modelled with the PEGASE code (Fioc & Rocca-Volmerange ll997t . 
These SSPs are created in a star-burst lasting for 25 Myr. We use a Scalo IMF (Scalo 1986). Three different metalicities are 
shown in the plots: solid: Z = 0.008, long-dashed: Z = 0.02, short-dashed: Z = 0.04. 



The stellar mass of this population is normalised to 10 
solar masses. We use the PEGASE code to produce flux ta- 
bles of broadband colours at linear spaced time-steps of 25 
Myr. This is done for the three different metalicities (Z = 
0.008; 0.02; 0.04), which are the input parameters into our 
code, where the chemical enrichment is approximated by linear 
interpolation of populations with these fixed metalicities. 

Fig Hshows the colour index evolution of our SSPs for the 
three different metalicities. 

The normalised luminosities in different broadband filters 
are shown in Fig.|4] and the mass-to-light ratios are shown in 
Figure [5] These SSP lookup-tables are used to assemble the 
total fluxes and stellar masses for a given star formation history. 
In Fig.|5]the fraction of initial mass in luminous stars, remnants 
and the sum of both, which goes into the stellar disc mass and 
the stellar mass-to-light ratio, is plotted for solar metalicity. 

4. Radiative transfer and the fitting procedure 

Here we describe the computation of the surface density pro- 
files and explain the main effects on the vertical profiles due to 
the variation of the disc parameters. 



1 1 


i i 


1 1 1 
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- N>... 








u 
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1 1 1 



2 4 6 8 10 12 14 



age [Gyr] 

Fig. 4. Luminosities of single stellar populations modelled 
with the PEGASE code with a Scalo IMF (Scalo 4986$ for solar 
metalicities. The luminosities are normalised to 1 solar mass. 
The different lines show different broadband filters: solid: U, 
long-dashed: B, short-dashed: V, dotted: R, dot-dashed: I. 
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Fig. 5. Evolution of the stellar mass-to-light ratios (including 
remnants) and of the mass of single stellar populations mod- 
elled with the PEGASE code with a Scalo IMF (Scalo [T98^ . 
The upper plot shows the M/L evolutions of a solar metalic- 
ity SSP in different broadband filters: solid: U, long-dashed: B, 
short-dashed: V, dotted: R, dot-dashed: I. We use the solar lu- 
minosities U = 5.54, B = 5.38, V = 4.75, R = 4.22, 1 = 3.87. The 
lower plot shows the mass evolution of the stellar population, 
normalised to its initial mass. The solid line shows the total 
stellar mass, and the other lines its two components separately: 
long-dashed: luminous stars, short-dashed: stellar remnants. 

The disc model is derived by performing a three-step it- 
eration cycle until good agreement of the model profiles with 
the observed surface-brightness and colour index profiles is 
reached. In the first step the normalised vertical disc profile 
is derived by adopting a star formation history and dynamical 
heating function combined with the relative contribution of gas 
and DM halo component. The intrinsic vertical profile is con- 
structed to fit the galaxy properties at the scaling radius i?o- In 
the second step the model is scaled with central surface den- 
sity, radial and vertical scale length of the stellar disc. Then the 
inclination with respect to the sky and the dust parameter are 
chosen to calculate the surface brightness distribution of verti- 
cal cuts parallel to the minor axis with a radiative transfer code. 
The geometrical parameters are varied until a good match to the 
set of V-band profiles is reached. Here the scaling of the gas and 



halo properties are also corrected iteratively. In the third step a 
multi-colour analysis is performed to analyse the deviations of 
the model to the complete set of observed surface-brightness 
and colour index profiles. These systematic deviations are used 
to estimate corrections for the star formation history and the 
heating function for an improved model. In the following these 
three steps are described in more detail. 

4. 1. The normalised vertical profiles 

We start with a pair SFR(t a — t) and a(t) expected for a Sc- 
galaxy from a basis set of normalised template functions of dif- 
ferent types. For the calculation of the self-consistent profile we 
need additionally the parameters of the gas and halo component 
normalised to the stellar parameter. These are 

1) the fractional surface densities of gas Q g and DM halo Qh 
(leading to Q s = 1 — Q g — Qh for the stars) 

2) the gas scale height normalised to the scale height of the 
stars Zd/z s determined implicitly by the maximum 'age' of the 
gas t g 

3) the velocity dispersion of the DM halo normalised to the 
maximum velocity dispersion of the stars «h = Ch/ce- 

Since Q g , Q^, Sh., and t g depend on the final scaling of 
the stellar disc to the observations, they must be corrected it- 
eratively. To fix the vertical profile we use the disc properties 
at the scaling radius Ro = lOkpc. For the gas component the 
surface density is derived from the total gas mass using the 
adopted radial scale length. The maximum 'age' t g is chosen 
to match the vertical scale length of the dust z g — z,i. For the 
dark matter halo we use the mass model of Sackett et al. ( 1994 1 
with disc M/L = 1 in a spherical DM halo. From the circular 
velocity of the halo u c h = 160km/s and with = v% h the 
halo velocity dispersion is = llOkm/s and the halo density 
is p h0 = vlJ(4nGRl) = 4.7 x 10" 3 M Q pc" 3 leading to 
the surface density Eh = 2z max /?ho = 47.5 Mqpc -2 up to a 
maximum height of z max = 5.1 kpc. The initial values for the 
stellar disc are also taken from the total stellar mass in Sackett 
et al. ( 1994 1 and then iterated. The normalisation of the veloc- 
ity dispersions are given by a c also determined iteratively from 
Eq.d] 

The calculation of the self-consistent profile includes the 
intrinsic luminosity profiles in all bands using stellar popula- 
tion synthesis. 

4.2. Extinction and projection to the sky 

The next step is to choose the stellar scale height z s (Eq. 1 1 81 . 
the radial scale lengths R s and the cutoff radius R max (Ea. l23> . 
the central surface density of stars So (Eas. ll6llT7l and the dust 
properties (Eq. l24> including the inclination i of the disc. 

With these parameters the projection of the galaxy to the 
sky is calculated. We use a simple radiative transfer code 
neglecting scattered light. Some test runs performed by M. 
Xilouris has shown that scattered light contributes only a few 
percent to the surface brightness of typical edge-on galaxies. 
This is different to the significant influence of scattered light on 
the determination of the dust temperature in the disc (Misiriotis 
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et al. 1200 1> . With the a;-axis along the line of sight, y-axis 
along the major axis, and z-axis along the minor axis we in- 
tegrate for a set of vertical cuts highly resolved z -profiles of 
surface brightness profiles in U, B, V, R, and I. The resolution 
is dxs'm(i) = 125 pc along the line of sight and dzsin(i) = 
12 pc parallel to the minor axis and we restrict the integration to 
the cylinder bounded by i? max and z max . Possible foreground 
extinction due to dust outside i? max is neglected. 

The extinction and reddening features of the vertical pro- 
files are very sensitive to the dust geometry and the inclination 
(for a detailed discussion see Just et al. Q596). Therefore we 
test a large range of parameters to find the best match to the 
observed profiles. 

Even if the match of the V-band profiles is not satisfactory 
we proceed to the multi-colour analysis to determine how to 
change the star formation history and/or the heating function. 

4.3. Multi-colour analysis 

Here we give an overview over the main effects on the vertical 
profiles due to the variation of the disc parameters. For details 
see Just et al. ( 1996). The global parameters of the stellar disc 
are determined by the surface brightness distribution at large 
heights z, where extinction can be neglected. The vertical cuts 
near the midplane are influenced by the dust distribution and 
the detailed distribution of the stellar subpopulations of differ- 
ent age. 

We first discuss the influence of the dust distribution. If the 
dust scale height is larger than the scale height of the intrinsic 
star light near the midplane, then extinction leads to a flatten- 
ing of the profile in the centre but no minimum (dust lane). 
Increasing the inclination results in a smearing out of the dust 
feature. Increasing the dust mass leads to a deeper and wider 
dust lane in the surface-brightness profiles but only a broaden- 
ing in the colour index index profiles; also the radial extension 
becomes larger. With a radial scale length larger than the stellar 
scale length the extinction feature is distributed more shallowly 
in the radial direction and the dust lane shows a larger offset to 
the major axis. 

The systematic variations of the shape of the vertical pro- 
files from red to blue are not only a signature of reddening near 
the midplane, but are also due to the relative distribution of 
young and old stars. Changing the star formation history leads 
to a different fraction of luminous young blue subpopulations. 
A higher fraction of young and intermediate age populations 
with 0,B and A-stars gives an excess in the surface-brightness 
profiles near the midplane which cannot be suppressed fully 
by dust extinction (especially in the red colours). The heating 
function determines the vertical spread of the subpopulations. 
Fast heating for the young populations results in centrally flat- 
tened surface brightness profiles combined with an extended 
blue regime. 

As an example of the differences in the intrinsic luminosity 
and colour profiles we show in Fig. ^] a comparison model 
with constant SFR (large fraction of young stars) and the 
heating function of the solar neighbourhood (quick heating of 
young stars). The effect is a reduced central luminosity and a 



Table 1. The observational data of NGC 5907. Most columns 
are self explaining. FWHM is the full width half maximum 
of the seeing. The two last columns give the total magnitudes 
and errors from our measurements, corrected for galactic fore- 
ground extinction. 



Filt 


Date 


Exp.time 


Airmass 


FWHM 


total 


phot 




Jul97 


[min] 






[arcsec] 


Magn 


error 


U 


11 


8 x 20 


1.10 - 


1.70 


1.5 


11.567 


±0.070 


B 


6 


8 x 20 


1.12 - 


1.53 


1.5 


11.279 


±0.054 


V 


5 


4 x 20 


1.06- 


1.13 


1.5 


10.419 


±0.057 


R 


5 


20 + 3 x 10 


1.30 - 


1.50 


1.3 


9.701 


±0.064 


I 


5 


4 x 10 


1.17- 


1.26 


1.8 


8.763 


±0.084 



very extended blue regime. Therefore less dust is needed in the 
central region to obtain the observed surface brightness pro- 
files. But the blue wings cannot be hidden by a reliable dust 
distribution. 

5. Observations and data reduction 

It is important to check the viability of our model by compari- 
son with multi-colour observational data of NGC 5709. In this 
section we describe the observations. The comparison will be 
presented in Sect. [6] 

The observations were performed during July 5-11, 1997 
using the 2.2m telescope of the Calar Alto Observatory, Spain. 
The CAFOS focal reducer was used changing the focal ratio of 
the telescope from f/8 to f/4.4. The CCD camera was equipped 
with a SITE 2048 x 2048 pixel chip. The pixel size was 24/Lt, 
corresponding to 0.53 arcsec. The circular field-of-use had a 
diameter of 16 arcmin. Standard Johnson U,B,V and Cousins 
R,I filters were used. Several images were exposed in each filter 
without interrupting the tracking of the telescope. The exposure 
times in UBVRI are listed in Tabled 

Sky flats were exposed in every dawn and dusk phase. 
The photometric quality was good or very good in all nights. 
Photometric standard stars (« 40) in the globular cluster M 92 
were exposed every night (Christian et al. 119851 Sandage & 
Walker lT966l Cafhev [T974l Stetson & Harris fT988b . 

5. 1. Data reduction 

The MIDAS program system was used for the data reduc- 
tion. From comparison and evaluation of all bias exposures an 
optimal working bias image was constructed. Similarly, from 
comparison and evaluation of all flat-field exposures optimal 
working-flat-field images were constructed for each colour. 
Bias subtraction and flat-field division were performed in the 
usual manner for each image. The flatness of the sky was 
checked in each image and corrected interactively, if necessary. 
The flatness of the sky down to less than 0.5% is an important 
condition for a reliable model of the surface-brightness distri- 
bution. 

The images in each colour were added together to a total 
image. A rebinning was not necessary since the tracking of 
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Fig. 6. Deep V image of NGC 5907 with the stellar mask. The image is rotated for a horizontal orientation of the galaxy. The field 
of view is 850 x 320 arcsec. The vertical flux profiles were calculated by mean values over the 13 indicated stripes of 45 arcsec 
width. The black lines mark the borders of these stripes. The flux values for the each point of the profiles were calculated by 
averaging along horizontal pixel rows of 45 arcsec length. The numbers in the stripes give their corresponding central positions 
in arcsec. NW is left, SE is right. 



the telescope was not interrupted between the exposures (this 
was checked, cf the column FWHM in Table [TJ. By stack- 
ing the images in the vertical direction it was possible to re- 
move the cosmic pixels during the summation (Kappa-sigma- 
summation). Fig.^shows as an example the total image in R. 
The faint filament extending northwards from the left part of 
the disk is real. It probably belongs to a remnant of an interact- 
ing dwarf galaxy (Shang et al. U996> . 



The flux calibration was calculated using the standard stars. 
The zero-points could be determined with a error between 0.05 
(B) to 0.08 (I) mag. The atmospheric extinction of the cor- 
responding airmass was calculated for each single exposure. 
Each exposure was considered with a corresponding weight 
in the flux calibration of the total image. The galactic fore- 
ground extinction for NGC 5907 is rather small: A g (B) = 0.01 
taken from RC3 (de Vaucouleurs et al. 1199 li t, see also in 
Burstein & Heiles ( 19841). The extinction corrections for the 
other filters were computed by the coefficients given in Cardelli 
et al. (1989 1 (their Table 3). The resulting corrections were 
A g = 0.012,0.010,0.008,0.006,0.004 mag for U, B, V, R, 
I, respectively. Because of the small distance of NGC 5907 the 
K-corrections are very small and were not considered. 



For the total magnitude of NGC 5907 we obtain the val- 
ues given in Table ^ These values are corrected for galactic 
foreground extinction. The errors are between 0.05 and 0.08 
mag. The corresponding B value from the RC3 is -Btot = 
11.11 ± 0.10. Our value is 0.17 mag weaker. The difference 
might be explained by the fact that the RC3 value was obtained 
by aperture photometry while we have masked the foreground 
stars in front of and around the galaxy. 



5.2. Surface-brightness and colour index profiles 

For the comparison with the models we need vertical profiles 
of the calibrated flux in UBVRI and the corresponding colour 
indices. All arithmetic calculations for this were done with the 
flux images. The flux calibrations were done afterwards. 

The total images were rotated to set the galaxy into a hor- 
izontal position. A mask with appropriate squares around all 
fairly bright foreground stars and background galaxies was 
constructed from the rotated total V image. 

For the construction of vertical flux profiles a num- 
ber of corresponding vertical stripes were defined. The cen- 
tre of the first stripe was placed at the centre of the 
galaxy and the others were placed at the horizontal positions 
±45, ±90, ±135, ±180, ±225, ±270 arcsec along the disk. 
(Fig. [6ji. This corresponds to 0.5, 1.0, 1.5, 2.0, 2.5, and 3.0 scale 
lengths. Each stripe had a width of 45 arcsec, thus the whole 
disk was covered. In each stripe we calculated the mean flux 
along each horizontal pixel row of 45 arcsec length. Whenever 
a stellar mask was hit, this was adequately considered in the 
mean value. In this way we obtained 13 flux profiles over z for 
each colour U, B, V, R, I. Further on we calculated the average 
of the left and right profiles of the corresponding horizontal 
positions. 

Using our flux-calibration formula these profile were trans- 
formed into magnitudes. The profiles of the colour indices U- 
B, B-V, V-R, R-I were obtained by subtraction of the corre- 
sponding profiles. 

In Fig.Qwe show the vertical cuts in U- and V-band as ex- 
amples. The thin lines are the NW- and SE-profile and the thick 
lines are the averaged profiles which we used for the model. 
The extinction feature (perturbed by the clumpy structure of 
the young subpopulation) is clearly visible. In the outer parts 
the profiles are exponential and the noise level is well defined. 
The small offsets of the profiles in the different profiles show 
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that the warp discussed in Morrison et al. (1994) seems to be 
not very systematic or is only significant for the young clumpy 
disc component. 

6. Results 

We show the comparison of the final model to the data in the 
surface-brightness and colour index profiles. Then we give the 
global parameters of the disc model, describe the scaling pro- 
cedure, and discuss the intrinsic properties of the model. 

6. 1. Comparison with the data 

In Fig. [9] the observed and model vertical surface-brightness 
profiles in all five bands are shown. The five profiles in each 
plot are vertical cuts parallel to the minor axis at distances 
along the major axes of 0, 2.4, 4.8, 7.2, 9.6, 12.0, and 14,4 kpc 
from top to bottom, respectively. The general features in the 
disc-dominated parts are reproduced reasonably well by the 
model in all bands. In the inner two profiles the additional 
bulge light can be observed and far above the galactic plane 
an additional more extended component appears as discussed 
in Lequeux et al. if 1996 1. Only the U-band profiles show a sys- 
tematic deviation in the sense that they seem to be radially more 
extended. The extinction features near the midplane are of a 
similar strength as observed but due to the strong clumping not 
reproduced in detail. The large optical depth totally obscures 
the strong excess stellar light in the intrinsic profiles (see Fig. 

EJi. 

The corresponding colour index profiles (Fig. 1101 are rel- 
atively flat as expected even at high optical depth due to the 
mixture of stellar emissivity and dust extinction along the line 
of sight (see Just et al. 1996 1. The systematic shape variation 
of the reddening features in the colour indices U-B, B-V, V- 
R, and R-I running from asymmetric with blue dip in the blue 
colours to more symmetric and redder in the red colours is also 
reproduced. 

6.2. Global properties 

The disc parameters of the final model are collected in Table|2] 
The inclination i = 87.5° and the cutoff radius i? max = 
19.7kpc are very similar to the model of Xilouris et al. J1999i . 
The radial and vertical scale length R s = 5.5 kpc and z s = 
520 pc are somewhat larger, because we concentrate our fit on 
the height regime \z\ < 1.5 kpc, where the Xilouris et al. fit is 
not satisfactory. 

For the total stellar mass (luminous plus remnants) we ob- 
tain from the central surface density S s ,o = 115Mqpc~ 2 a 
total stellar mass of 

M s = 1.9 x 10 10 Mq (27) 

corresponding to an average star formation rate of (SFR) = 
M s / <7eff Aa = 2.34 MQ/yr and a present day value of SFRq = 
1.64M /yr. 

The dust component has a radial and vertical scale length 
of i?d = 7.7 kpc and z^ = llOpc, respectively. The central ex- 
tinction coefficient is Ay.o = 4.0mag/kpc, which corresponds 



stellar disc dust component 



i 


87.5° 






flmax [kpc] 


19.7 






fl s [kpc] 


5.5 


R d [kpc] 


7.7 


Zs [pc] 


520 


Zd [pc] 


110 


So [Mq pc" 2 ] 


170 


A v , [mag/kpc] 


4.0 


ffoff 


0.676 






a 


1.05 







Table 2. Parameters of the final model for NGC 5907. At the 
left hand side the physical quantities of the stellar disc are given 
and at the right hand side the parameters for the dust compo- 
nent. Inclination i and cutoff radius i? max are the same for all 
components. So is the integrated star formation at the centre 
and g ff is the conversion factor to the stellar surface density 
^s,o = ffeff^o = 115 Mqpc -2 . C z is the shape correction fac- 
tor of the stellar density profile. 



to a face-on optical depth of t v = 0.81 and central gas den- 
sity of pd.o = 7.0 x 10 -4 Mq/pc 3 . The total dust mass is 
Md = 47ri? 2 ZdPd,o = 5.7 x 10 7 Mq leading to a gas-to-dust 
ratio of F = M g /M d = 209. 

6.3. The scaling radius R 

The input parameters of gas and halo component for the intrin- 
sic disc structure are determined at the scaling radius Rq = 
10 kpc (see Tabled For the normalisation we need the proper- 
ties of the stellar component from the final model. The stellar 
surface density is S s = 18.7 Mqpc~ 2 (Ea. l23t and the maxi- 
mum velocity dispersion is a e = 32km/s (Eq.GU. 

For the gas component we adopt the same scale lengths as 
for the dust, i.e. R g = i?d = 7.7kpc and z g = z d = llOpc, 
which for the radial distribution is consistent with the larger ex- 
tension of the HI distribution. The vertical scale height corre- 
sponds to the maximum age bin of t g = 3 Gyr. This can be con- 
verted to a maximum velocity dispersion of er g m = 0.325er c = 

10.4 km/s and an rms value of a g = 0. 16ct = 5.1 km/s for the 
gas. 

For the DM halo the surface density and the velocity disper- 
sion taken from the mass model of Sackett et al. ( 1994 1 Eh = 

48.5 Mqpc~ 2 up to a maximum height of z max = 5.1 kpc and 
o"h = 110 km/s (see Sect. 14. 1> are used. Since we deduced the 
halo properties from the rotation curve, we checked the recon- 
struction of the rotation curve with our final model. A mass 
model leading to an equally good fit of the rotation curve in 
Sackett et al. can be constructed by a four component galaxy 
including a Hernquist bulge, exponential discs for the stellar 
and the gas component, and a pseudo-isothermal dark matter 
halo. The bulge with a mass of 3 x 10 10 Mq and a core ra- 
dius of 1 kpc reproduces the steep rise of the CO rotation curve 
presented in Sofue et al. {1996 1. We adopt a flattening of 0.5 
as determined in Barneby & Thronson ( 1994). The exponential 
discs are determined by the scalelengths of 5.5 kpc and 7.7 kpc 
and the total masses of 1.9 x 1O 1O M and 1.2 x 1O 1O M 
(corresponding to 2.18 x 10 10 Mq and 1.65 x 10 10 Mq when 
neglecting the cutoff) for stars and gas, respectively. The dark 
matter halo has a local density of 5.5 x 10 _3 Mq/pc 3 at 
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Fig. 7. Vertical surface-brightness profiles of NGC 5907 each averaged over 45" in the radial direction with an offset to the 
minor axis of 0, 45, 90, 135, 180, 225, and 270" from top to bottom. For a better visibility the magnitude scale is shifted by 2 
magnitudes from one to the next. Long and short dashed lines give the individual profiles for the NW(-) and SE(+) sides of the 
disc and the thick full lines are the averaged profiles. 




offset arcsec offset arcsec 

Fig. 8. Vertical profiles of the colour indices U-V and B-R for NGC 5907 for the same cuts as in Fig.0 From top to bottom, the 
magnitude scale is shifted in colour by -0.5 magnitudes from one to the next. 



i?o = lOkpc with a core radius of 1 kpc. The local density 
is slightly higher than that used in our model in order to cor- 
rect for the effect of the core. The rotation curve of NGC 5907 
presented in Sofue et al. {1996 ) is significantly different and 
also shows the uncertainties in the observational data itself. To 
reproduce the prominent maximum around a radius of 12 kpc, 
one can use a NFW halo with scale radius of 6 kpc and lo- 
cal density of 4.78 x 10 _3 Mq/pc 3 (as in our model) com- 
bined with a smaller bulge with core radius 0.5 kpc and mass 



of 1.2 x 10 Mq. There is a large freedom in the parameter 
choice to reproduce the rotation curves, especially if we allow 
the disk mass to vary. Instead of using a NFW profile the max- 
imum in the rotation curve of Sofue et al. can also be reached 
by increasing the disc mass by a factor of 2-3. 

The total surface density at i?o is Stot = 76.0 Mqpc~ 2 , 
which is used to determine the ratios Q s = £ s /£ to t etc. The 
halo velocity dispersion is normalised by Sh = Oh/oc = 3.5. 
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Fig. 9. Comparison between observations and model. Vertical surface brightness profiles for NGC 5907 in the U, B, V, R, and I 
bands from upper left to lower right. The seven profiles in each plot are vertical cuts at distances along the major axes of 0, 45, 
90, 135, 180, 225, and 270" corresponding to 0, 2.4, 4.8, 7.2, 9.6, 12.0, and 14.4 kpc from top to bottom. Dashed lines represent 
the data, while solid lines give the model. From top to bottom, the magnitude scale is shifted downwards in brightness by 0, 2, 4, 
6, 8, 10, and 12 magnitudes, respectively. 



at Ro — 10 kpc 



scaling parameter 





520 


<j c [km/s] 


32 


E s [M pc- 2 ] 


18.7 


Qs 


0.25 




8.8 


Q s 


0.11 


s h 


48.5 




0.64 


o- h [km/s] 


110 


Sh 


3.5 


2g[pc] 


110 


i g [Gyr] 


3.0 



Table 3. Input parameters for the final model of NGC 5907. At 
the left hand side the physical quantities at Rq = 10 kpc are 
given and at the right hand side the corresponding normalised 
parameters of the normalised profile. The corresponding de- 
scription is given in Sect. 16. 31 



All input parameters of the final model for NGC 5907 are col- 
lected in Table [3] 



6.4. Intrinsic disc properties 

The self-consistent vertical profiles of the final model are 
shown in Fig. The stellar density profile is approximately 



exponential with an excess below \z\ « 500pc. At the mid- 
plane the stellar density compared to the extrapolated expo- 
nential profile shows an excess of a factor of 2.1. The cor- 
responding excess in surface density is 48%. The gas den- 
sity profile matches the scale height of the dust component 
with z g — Zd — llOpc. The central gas density is slightly 
(5%) larger than the central stellar density. The DM halo pro- 
file is very flat and decreases by 44% from the midplane to 
Zmax = 5.1 kpc. The halo density exceeds the disc density 
(stars plus gas) for \z\ > 510 pc, but the cumulative mass, 
which determines the local gravitational force, dominates only 
at \z\ > 2.7kpc. 

The left panels in Fig.El s h° w the intrinsic luminosity pro- 
files in all bands. For comparison the density profile is also 
plotted to emphasise the enhanced stellar light of the young 
stars near the midplane. As the lower panel shows for the B- 
band, the profiles deviate from the exponential shape already at 
heights below \z\ = 1 kpc. In the right panels the corresponding 
colour profiles and mass-to-light ratios are shown. The intrin- 
sic colours above \z\ ~ 200 pc show very shallow gradients for 
this model. The colour gradient due to the blue young popula- 
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Fig. 10. Vertical profiles for NGC 5907 in the U-B, B-V, V-R, and R-I colour indices from upper left to lower right. The seven 
profiles in each plot are the same vertical cuts as in Fig.|9] Dashed lines represent the data, while solid lines give the model. From 
top to bottom, the magnitude scale is shifted by 0, 1, 2, 3, 4, 5 and 6 magnitudes to the blue, respectively. 
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Fig. 11. Vertical density profiles of the final of NGC 5907 
model normalised to the values at the scaling radius Rq = 
lOkpc. The gas density (short dashed line) slightly dominates 
towards the midplane. The stellar density (full line) has an 
exponential scale height of z s — 520 pc with a density en- 
hancement to the midplane. The halo density (long dashed line) 
shows only 20% variation below 2.5 kpc due to the adiabatic 
contraction. 



tion can be seen only in the innermost part. The mass-to-light 
ratio is more sensitive to the age-distribution of the stars and 
varies much more strongly than the colours. It varies by a factor 
of 10 in U-band and even in I-band by a factor of 3. The M/L- 
gradient is significantly larger than zero up to \z\ ~ 500 pc at 
least in the blue bands. The average value of M/L corresponds 
to the local value near a height of \z\ ~ 150 pc (dots in the 
lower right panel of Fig.ll2i. 

The luminosity weighted mean metalicity of the final model 
varies from [Fe/H] = 0.06 at the midplane to [Fe/H] = 
—0.59 at high latitude \z\, which is only slightly less than the 
variation in the chemical enrichment of the subpopulations. 
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(M/L) cxp 
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2.39 
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Table 4. Intrinsic properties of the final model of NGC 5907 
in the different bands. The absolute quantities (central emis- 
sivity, surface brightness and total luminosity) depend on the 
central stellar surface density and the radial and vertical scale 
length as given in Table|2] The first row gives the central emis- 
sivity eq and the second row shows the enhancement factor 
(e/e eX p)o with respect to the exponential extrapolation form 
z w 3z s to the midplane. The mass-to-light ratios (M/L)q at 
the midplane are very low due to the high fraction young stars. 
The next two rows give the central face-on surface brightness 
L(R = 0) in physical and observational units, mtot is the total 
apparent luminosity with the distance module Am = 30.2 mag 
and Am tc = m to t — "icxp quantifies the additional stellar 
light compared to pure exponential vertical profiles. The aver- 
age stellar mass-to-light ratios (M/L) to t of the model are the 
same for each vertical profile and the whole disc. For the ex- 



ponential extrapolation (M/L) cxp is constant all over the disc 
and equals the MJ L of the model at large heights z. 



In order to show that the evolutionary history cannot be de- 
composed arbitrarily into a dynamical heating function and an 
appropriate star formation history, we calculated comparison 
models, where we have fixed the heating function by that of the 
solar neighbourhood (Fig. |2]lower panel, thin line) and changed 
the SFR in order to get a higher fraction of young stars. But 
even with a constant SFR the strong heating of the young sub- 
populations spread the blue stars to a much larger height yield- 
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Fig. 12. Upper left panel: The emissivity of the central profile of NGC 5907 in magpc~ 3 . For comparison the stellar density 
profile is also plotted in units of Mq pc~ 3 /10 (also in mag). Lower left panel: The B and I-band emissivity for comparison 
with a model using the heating function of the solar neighbourhood and a constant star formation rate. For the B-band profile the 
exponential profile fit (dot-dashed line) is shown. The corresponding density profiles are normalised to Mq pc _3 /5. Upper right 
panel: Intrinsic colour profiles at the centre of NGC 5907 (thick lines). The profiles are shifted to avoid overlap. The profiles of 
the comparison model (thin lines) with constant SFR show more extended colour gradients due to the faster dynamical heating 
(see Fig. |2j. Lower right panel: The mass-to-light ratios in all bands, also shifted to avoid overlap. The dots denote the average 
values for the total surface density (see Table0}. 



ing more extended colour gradients and a smaller stellar light 
excess near the midplane (thin lines in Fig.ll2i. 

7. Discussion 

7. 1. Disc evolution 

Our final model combines a star formation history with a mod- 
erate maximum at an age of t — 7.2 Gyr with a slow dynamical 
heating for the young stellar population compared to the so- 
lar neighbourhood (see Fig.|2li. If we use as in the comparison 
model the heating function of the solar neighbourhood (con- 
stant increase of kinetic energy of the subpopulations with age), 
a higher present day SFR is necessary to get a comparable lu- 
minosity in the blue bands near the midplane. Even a constant 
SFR (which is often claimed even for the Milky Way) does 
not reach the intrinsic luminosity of our final model, but leads 
to additional blue light at intermediate heights \z\, where dust 
reddening cannot compensate. We do not claim that the cho- 



sen SFR is unique, because the parameter space is too large 
and a model on a much higher level of comparison including 
radial variation of the intrinsic profiles and clumping near the 
midplane would be necessary. Our result shows that a conti- 
nously decreasing SFR in the last few Gyr can account for 
the observed high intrinsic luminosity at the midplane neces- 
sary for the strong NIR-submm brightness. A recent star burst 
is not needed but also not excluded by our model at this level 
of modelling using a simple radial disc structure. 

7.2. Mass-to-light ratio 

The stellar mass-to-light ratio Mj L strongly depends on the 
IMF at the low mass end without changing the luminosities 
much. In our model with the Scalo-IMF we get a relatively low 
value M/Ly = 0.91, which fits well to the mass models of 
Sackett et al. ( 1994). Since the disc of an Sc galaxy is younger 
than an Sbc galaxy like the Milky Way, this value also fits well 
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compared to the values in the solar neighbourhood. A rough 
estimate of the excess in blue light of the young stellar popula- 
tion in our model compares well to the additional light needed 
in Misiriotis et al. (2001 1 to heat the dust. Their corresponding 
mass-to-light ratio and recent SFR are slightly different, be- 
cause they used a different stellar population model (Salpeter 
IMF and no mass loss due to stellar evolution). 

The low mass-to-light ratio of the stellar disc is a clear in- 
dication of a low disc mass model as discussed also in Sackett 
et al. d!994i . Our underlying mass model is consistent with 
the observed rotation curve of Sackett et al. with an Pseudo- 
isothermal DM halo and for the rotation curve of Sofue et al. 
( 19961 a NFW halo fits better. But since the mass model is fully 
self-consistent at the scale radius i?o = lOkpc only, we cannot 
distinguish between isothermal or NWF profiles for the dark 
matter halo. 

7.3. Dust component 

For the dust component we find a similar scale height as 
Misiriotis et al. (2001 1 and a similar total dust mass, which is 
responsible for the FIR brightness. A significant difference to 
the models of Xilouris et al. dl999> and Misiriotis et al. is the 
larger dust scale length of i?d = 7.7 kpc. The effect of a larger 
i?d is, apart from the stronger extinction feature at large dis- 
tance form the minor axis, a shift of the extinction dip to higher 
vertical offset at all vertical cuts. This effect cannot be substi- 
tuted by a larger inclination i, because this would smear out the 
extinction feature to larger heights above the major axis. A side 
effect of the larger radial scale length is the reduced face-on op- 
tical depth at the centre (Yy = 0.81 instead of 1.5). Xilouris et 
al. also found that the dust scale length is larger than the stellar 
scale length in their best fit. The large radial scale length of the 
dust also fits to the more extended distribution of HI and of the 
FIR brightness. 

8. Conclusions 

We have observed deep photometry of NGC 5907 in the U, B, 
V, R, and I-band in order to construct a physical disc model 
based on a multi-colour analysis. We have constructed a self- 
consistent evolutionary disc model for NGC 5907, which fits 
the general characteristics of the surface brightness distribu- 
tion in all bands. The stellar disc is built up by a sequence of 
stellar subpopulations distributed according to a smooth SFR 
and dynamical heating function a(t). The emissivity is calcu- 
lated with stellar population synthesis and the appearance of 
the galaxy is calculated by radiative transfer including an ex- 
ponential disc component of dust. 

This is the first physical model of the disc of NGC 5907 that 
gives the correct luminosities and colours of the outer exponen- 
tial parts of the disc and also the large amount of dust observed 
in FIR/submm. The reason for the corresponding excess in stel- 
lar light near the midplane of the disc is the concentrated lumi- 
nous young subpopulation produced by the continuous SFR. 
Our model is a significant improvement of the disc model of 
Misiriotis et al. (2001 1, who combined the phenomenological 
exponential model of Xilouris et al. 0999 1 with an additional 



thin disc of uniformly mixed young stars and dust in order to 
produce the excess in FIR/submm luminosity. 

This evolutionary disc model naturally explains three dif- 
ferent aspects known from the literature: 1) The exponential 
vertical structure well above the midplane as modelled phe- 
nomenologically by Xilouris et al.( 1999 1 is the result of the 
older subpopulations in dynamical equilibrium in the total 
disc potential. 2) The large amount of dust observed in the 
FIR/submm (modelled by Misiriotis et al. 2001) is necessary 
to hide the bright young subpopulations confined to the mid- 
plane. The young stars are the smooth end of the SFR without 
a recent star burst. 3) The additional light necessary to heat the 
dust is automatically delivered by the high mass-to-light ratio 
near the midplane. 
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